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ABSTRACT 


In  is  report,  we  discuss  accurate  and  robust  sliding  mode  tracking  control  for 
highly  nonlinear  robot  manipulators  using  a  disturbance  observer.  To  eliminate  the 
chattering  problem  existed  in  conventional  Sliding  Mode  Control  (SMC)  approach, 
which  is  caused  by  modeling  errors  and  uncertainties.  The  efficient  compensation 
of  the  disturbance  observer  has  been  introduced. 

The  proposed  sliding  mode  control  is  presented  in  two  theorems.  The  bounded 
stability  of  the  proposed  control  method  is  proven  and  the  efficiency  of  the  control 
algorithms  has  been  demonstrated  by  simulations  for  a  position  tracking  control 
of  a  two-link  robot  subject  to  parameter  and  payload  uncertainties. 


1  Introduction 


The  overall  robot  control  problem  can  be  divided  into  two  phases,  trajectory  planing  and 
control  phases.  The  trajectory  plcining  is  the  placing  of  the  desired  time  trajectory  from 
an  initial  position  to  a  desired  position  of  robot  manipulators  with  collision  avoidance  if 
necessary.  The  control  phase  is  the  tracking  control  of  robot  manipulators  to  track  the  de¬ 
sired  trajectory  planned  in  the  trajectory  planning  phase.  This  tracking  control  phase  is 
the  field  of  our  concern  in  this  paper.  The  strategies  for  designing  a  robust  and  accurate 
tracking  controller  for  a  highly  nonlinecir  robot  manipulator  have  been  studied  enthusiasti¬ 
cally  in  order  to  extend  its  application  fields.  There  are  several  approaches  to  attempt  to 
obtain  the  desired  tracking  performance  such  as,  Computed  Torque  Method  (CTM)  [1-3], 
.Adaptive  Control  [4,  5],  Sliding  Mode  Control  (SMC)  [6-13],  and  so  on  [14].  Each  method 
has  its  merits  and  shortcomings.  The  CTM  normally  provides  a  feasible  controller  if  tlie 
exact  knowledge  of  the  manipulator  dynamics  is  available.  However,  for  a  large  amount 
of  applications,  it  is  impossible  to  obtain  the  complete  dynamic  model  of  robots,  due  to 
modeling  uncertainties,  parameter  variation  and  unknown  payloads.  These  uncertainties, 
especially  the  error  of  inertia  matrix,  may  result  in  the  instability  of  robot  systems.  .All  of 
highly  nonlinear  model  dynamics  are  taken  into  account  in  order  to  calculate  the  control 
input  which  is  a  hard  load  in  view  of  the  computation  time.  Thus  the  optimal  trade-off  is 
made  between  approximate  modeling  of  the  robot  manipulator  and  its  output  performance. 
The  CTM  is  sensitive  to  modeling  errors  which  come  from  linearization  and  approximation 
o'"  nonlinearities,  uncertainties  in  physical  systems,  and  disturbances  [3]. 

In  order  to  give  the  robustness  to  modeling  uncertainties  and  parameter  variations,  adap¬ 
tive  algorithms  for  a  highly  nonlinear  robot  manipulator  eire  considered  such  as  Self-Tuning 
and  Model- Reference  Methods  but  essentially  nccd  an  adaptation  mechanism  for  the  pa¬ 
rameter  identification  resulting  in  a  heavy  computation  burden,  complexities,  and  a  high 
cost  for  the  digital  implementation.  As  another  robust  controller  different  from  the  adaptive 
approach,  SMC  using  Variable  Structure  Systems  (VSS)  for  robot  arms  is  studied  by  many 
researchers  recently  [6-13].  The  desired  performance  with  a  simple  control  structure  can 
be  obtained  in  the  existence  of  the  acceptable  model  error  and  unknown  payload  using  the 
"sliding  mode".  In  the  sliding  mode,  it  is  well  known  that  the  controlled  system  is  robust 
to  the  bounded  parameter  variations  and  disturbaxices  [17,  18].  The  other  advantages  of  a 
SMC  are  that  the  output  performance  can  be  predetermined  by  a  sliding  surface  and  there 
is  no  overshoot  in  regulations.  The  first  application  of  SMC  to  robot  manipulators  seems  to 
be  in  the  work  of  Young  dealing  with  a  set  point  regulation  problem  [6].  A  modification  of 
the  Young  controller  Wcis  presented  by  Morgan  [8].  Other  SMC  of  robot  manipulator  may  be 
found  [10-12].  The  SMC  unfortunately  has  a  problem  that  the  motion  trajectory  of  system 
changes  frequently  in  the  vicinity  of  the  sliding  surface  due  to  a  high  frequency  switching 
of  the  input  which  is  called  cis  "chattering" .  Because  in  order  to  guarantee  the  existence  of 
the  sliding  mode,  the  control  has  imperatively  the  discontinuities  designed  by  the  maximum 
bound  of  model  errors  or  uncertainties.  The  conventional  SMCs  for  robot  manipulators  use 
the  feedforward  compensation  by  using  an  available  model  dynamics,  there  exist  inevitable 
chattering  due  to  model  errors  or  uncertainties  in  system  parameters.  To  alleviate  this 
problem  bounded  layer  method  or  saturation  function  [9]  are  used  instead  of  discontinuous 
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parts.  These  methods  are  based  on  the  try-and-error  not  a  systematic  approach  and  not 
enough  to  obtain  the  desired  performance  with  the  satisfactory  tracking  error.  Thus  as  pos¬ 
sible  as  the  exact  compensation  of  the  high  nonlinear  interactions  in  robot  link  systems  is 
essential  to  obtain  the  good  performance  without  chattering.  Using  the  multirate  sampling 
concept,  a  SMC  with  am  off-line  feedforward  compensation  for  the  nonlinear  interaction  of 
robot  manipulators  using  an  available  model  is  presented. 

In  this  paper,  design  of  a  sliding  mode  control  algorithm  with  an  efficient  on-line  com¬ 
pensation  for  robot  manipulators  is  studied  for  robust  and  accurate  tracking  of  the  desired 
trajectories.  The  proposed  control  consists  of  two  parts  of  the  efficient  compensation  and 
the  sliding  mode  only  based  on  the  nominal  inertia  matrix  of  robot  manipulators.  This 
compensation  technique  is  responsible  for  not  only  highly  nonlinear  interactions  but  also 
modeling  errors  of  the  inertia  matrix  using  the  acceleration  information  and  the  nominal 
inertia  matrix.  The  sliding  mode  control  is  applied  to  the  compensated  dynamics  of  robot 
manipulators  to  obtain  robust  and  accurate  tracking  performance  with  the  correction  of 
compensation  errors  in  view  of  the  real  implementation  using  the  microprocessor.  The  sta¬ 
bility  of  the  proposed  SMC  is  deeply  analyzed  in  sense  of  Lyapunov,  summarized  in  two 
main  theorems  with  a  preliminary  lemma.  Based  on  this  analysis,  the  given  desired  track¬ 
ing  specifications  can  be  satisfied  by  a  proposed  SMC.  The  design  of  the  proposed  SMC  is 
independent  of  the  uncertainty  parameter  space.  Thus  there  can  be  no  chattering  due  to 
modeling  errors  or  uncertainties. 

The  rest  of  the  paper  is  organized  as  follows.  Section  2  gives  the  mathematical  basis  for 
a  stable  SMC  design  by  Lemma  1.  In  Section  2,  a  new  SMC  with  the  efficient  compensation 
is  proposed,  and  its  stability  is  analyzed  in  Theorem  1.  A  modified  control  is  then  presented 
in  Section  3  and  cinother  more  efficient  approach  is  proposed  in  Section  4.  Its  stability  is 
summarized  in  Theorem  2.  The  simulation  studies  for  the  proposed  algorithms  are  carried 
out  in  Section  5  and  conclusions  are  presented  in  Section  6. 


2  Mathematical  Work 

The  motion  equations  of  an  n  degree-of-freedom  manipulator  can  be  derived  using  the 
Lagrange- Euler  formulation  as 

‘^(<7(0,  <A)  •  qit)  +  D{q{t),q{t),(f>)  =  r(0  (1) 

where  ./(g(t),<^)  G  fT”'"  is  a  symmetric  positive  definite  inertia  matrix.  D{q{t),q{t),0)  G 
is  called  a  smooth  generalized  disturbances  (SGD)  vector  including  the  centrifugal  and  Corio¬ 
lis  terms  H {q{t) , q{t) ,  (f>)  G  W',  Coulomb  and  viscous  or  any  other  frictions  F(q{t),q{t).  o)  G 
/?",  gravity  terms  G{q{t),4))  G  R". 

D{q{t),q{t),(l>)  =  H{qit),q{t),(f>)  +  F{q{t),q{t),4>)  +  G{q{t),(i>) 

where  r(7)  G  R  is  an  input  vector,  q{t),  q(t),  and  q(t)  G  R  are  the  generalized  position, 
velocity,  and  acceleration  vector,  respectively,  (b  is  a  vector  composed  of  the  parameters  of 
robot  manipulators  (i.e.  the  masses,  lengths,  offset  angles,  and  inertia  of  links).  Exact  mod¬ 
eling  of  robot  dynamics  is  difficult  because  of  parameter  uncertainties,  unknown  frictions. 
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payloads  variations.  Therefore,  by  using  estimated  parameters  in  the  model,  we  express  the 
model  of  the  robot  system  (1)  as  follows 

q{t)  + =  r{t).  (2) 

The  objective  of  robot  trajectory  control  is  to  follow  the  desired  trajectories  qd(t),  qd{i)- 
qd{t)  computed  in  planning  phase.  Let’s  define  the  state  vector  X{t)  6  for  SMC  as 

X(0  =  [^i(0^  (3) 


where  Xi{t)  and  X2{t)  €  are  the  trajectory  errors  and  its  derivative 

^i(0  =  e(0  =  (qdit)  -  qit)) 


■  0  7  ■ 
0  0 

■X(t)- 

0 

•r(<)  + 

0 

J{q{t),<f>)-'^  ■  D{q{t),q{t),<f>) 

+ 

0 

.  . 

X2it)  =  e{t)  =  (qdit)  -  q{t)). 

Then  the  state  equation  of  the  robot  system  becomes 

.Y(f)  = 

An  augmented  sliding  surface  vector  s{t)  is  defined  as 

s{t)  =  X2it)  +  K,  ■  Xx{t)  +  K,-  J  X^{t)dt 

where  and  Kp  €  are  gain  matrices.  The  ideal  sliding  dynamics  defined  by 
(5)  is  the  trajectories  q',  q^,  q*  €  which  satisfy  the  following  equation 

x;{t)  +  K,  •  x;{t)  +  Kp  ■  x:{t)  =  0. 

.Y*(f)  =  A  •  X‘{t) 


(o) 

(0) 

Equation 


Or 

where 


(S) 


x-{t)  = 


{qd{t)-q'{t)V 
{qd{t)  -  qAt)VV 


€  /?•"*,  A  = 


0  I 
-Kp  -K, 


e  R 


2n  X  2n 


We  can  choose  and  Kp  so  that  all  the  eigenvalues  of  A  have  negative  real  parts,  which 
guarantees  the  exponential  stability  of  the  system  (7),  i.e.,  there  exist  positive  constants  I\ 
and  K  such  that 

(')) 

where  ||  ■  ||  is  the  induced  Euclidean  norm  [15). 

Now,  we  will  state  the  following  lemma  as  a  prerequisite  to  the  two  main  theorems. 
Lemma  .  :  If  the  sliding  surface  defined  by  Equation  (6)  satisfies  ||  s(t)  ||<  7  for  any 
t  >  to  and  II  X(to)  ||<  7/K:  is  satisfied  at  the  initial  time,  then 


(ina) 
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(10b) 


11^2(0  i|<£2 

is  satisfied  for  all  t  >  to  where  £i  and  are  positive  constants  defined  as  following 

£i  =  --7,  £2  =  7-[1  +  ^--],  Z=\\[K,K,]\\.  (10c) 

/C  /c 

Proof  [16] 

Define  a  new  state,  Xolt),  to  be  aji  integral  of  the  state  Xi{t),  i.e., 

>^o(0  —  /  Xi{u)du,  Xo(to)  =  0.  (11) 

Jto 

then  Xo(t),  the  augmented  sliding  surface  (6)  can  be  rewritten  as 

s(t)  =  ^2(0  "h  '  -^i(0  "h  '  Xo{t).  ( 12) 

If  we  define  a  new  state  vector  X{t)  =  [Xo(t)^,  Xi(t)^]^,  thus  we  have  the  following  state- 
space  equation 

X(0  =  A-X(0+  J  -sit).  (13) 

Where  s{t)  may  be  considered  as  bounded  disturbances,  [j  s{t)  1|<  7.  The  solution  of 
Equation  (13)  is 

J  •s(u)]du.  (M) 

From  Equation  (9)  and  boundedness  of  s(t),  the  norm  of  X{t)  becomes 

II  m\\<K-  II  X(t,)  II  +  f  II  II  7 

Jto  t 

II  XW  ||<:^.7  +  1||  X(Q  ||-2).A-.e-“"-' 

K  K 

<  7  •  —  for  all  <  >  to,  (15) 

K 

which  drives  Xi(t)  to  satisfy  the  following  inequality 

II  ^.(0  ll<-  T 

K 

From  Equation  (12)  we  have 

X2it)  =  3{t)^[KoKo]-Xit).  (16) 

then  by  using  the  result  of  Equation  (15) 

II  x,{t)  ||<  7  .  (1  +  Z  .  ^i| 

K 


_  gA(«-to) 


Xito)  + 
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for  all  t  >  to  which  proves  (10b)  and  therefore  we  complete  our  proof. 

Lemma  1  implies  that  the  tracking  error  of  the  trajectory  and  its  derivative  are  uniformly 
bounded,  provided  that  the  sliding  surface  is  bounded  for  all  time  t  >  to-  Using  this  results 
of  Lemma  1,  we  can  give  the  specifications  of  tracking  errors  dependent  on  convergence  rate 
which  is  determined  by  the  sliding  surface,  (9).  In  the  next  section,  we  are  going  to  design 
a  controller  with  efficient  compensation  which  guarantees  the  boundedness  of  s(<),  i.e.,  if 
II  •s(t)  11  <  7  for  a  given  7,  then  the  trajectory  error  is  bounded  in  virtue  of  Lemma  1. 


3  Efficient  Sliding  Mode  Control  for  Robot  Manipu¬ 
lators 

There  are  a  variety  of  disturbances  acting  on  robot  manipulators,  which  make  robust  tracking 
control  highly  desirable.  It  is  commonly  noted  that  the  generalized  nonlinear  disturbances. 
D{q{t),q{t),<f>),  should  be  compensated  to  improve  performance.  For  sliding  mode  control  of 
the  robot  system  (4),  the  equivalent  control  of  the  augmented  sliding  surface  (5)  is  obtained 
[17,  18] 

T^git)  =  D{q{t),q{t),<f>)  +  J{q{t),<l>)  ■  {qd{t)  +  +  K^Xi)  (17) 

The  SGD  D{q{t),q{t),<f>)  is  included  in  an  equivalent  control,  Te,(t).  SGD  is  generally  com¬ 
plex,  and  difficult  to  be  calculated  directly  due  to  uncertainties  in  the  model. 

In  the  paper,  using  the  efficient  estimation  method  so  called  disturbance  observer,  we 
consider  the  following  continuous  control  input,  r(<) 

T{t)  =  Tc{t)  +  T,{t)  (IS) 

Where  r<.{t)  is  the  compensation  term  for  SGD  as  well  as  the  error  of  nominal  inertia  matrix, 
is  not  the  direct  calculation  from  D{q{t),q{t),^)  in  the  model  (2)  but  the  efficient  estimation 
of  the  generalized  disturbance,  D{q{t),q{t),d>)  only  using  the  nominal  inertia  matrix.  ./%■. 
of  the  model  (2)  and  an  available  acceleration  information  which  can  be  obtained  from  the 
speed  information  by  means  of  Euler  method. 

Tcit)  =  r(t)  -  Jfsf  ■  q{t) 

=  D{q{t),q{t),(i>)  -t-  AJ{q{t),<i>)  ■  q{t)  -  J[q{t),<i>)  ■  Xq{t)  -  Ar(<)  (10) 

where  L\J{q{t),  (f>)  is  the  deviation  between  the  real  inertia  matrix  and  its  nominal  value. 
Xij  is  the  acceleration  inf  rmation  error  to  the  real  acceleration  value,  Xr(t)  is  the  control 
input  delay  error  resulting  from  the  digital  control  are  defined  by 

XJ{q({t),d>)  =  J{q(t),4>)  -  J/v 

Aq{t)  =  q{t)  -  q{t) 

Ar(f)  =  r(t  —  h)  —  r(t) 

where  h  is  the  sampling  time.  The  second  term  in  Equation  (16)  is  defined  as 

r,it)  =  (f^git)  +  r^it))  (20) 
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where  fgg{t)  is  the  modified  equivalent  control  for  the  compensated  dynamics  of  Equation 
(1),  and  is  designed  so  that  the  error  dynamics  of  the  controlled  system  has  the  sliding 
surface  dynamics  defined  by  Equation  (7).  It  is  defined  cis 

^eq{t)  ~  Jn  •  Kt,  ■  X2  +  Kp  ■  Xi).  (21) 

T^{t)  is  the  feedback  term  of  the  sliding  surface  resulting  from  the  compensation  error 

•5  f  ^  ) 

rx(<)  =  J/V  •  -^(0  +  ^x2  ■<7-a(0}.  ^“(0  = 

If  we  apply  the  input  control  torque  given  by  Equations  (18)-(22)  to  the  robot  system  (5). 
we  obtain  the  following  equation 

Xiit)  =  -J~\q{t),d>)  •  {AJ{qit),<f>)  ■  q{t)  -  Ar(t))  +  Aq{t)  +  q^ 

-  J~^{q{t),  0)  •  J/v  •  [qd  +  •  X2  +  Kp  ■  Xx  +  ky^x  ■  +  ^x2  '  7j — /  ,!  ,!  r]  (  -'5) 

11  II 

and  the  dynamics  of  s{t)  is  expressed  in  the  following  simple  form 

sit^ 

s{t)  =  n,(t)  -  [^^1  -  sit)  4-  k^2  ■  II  II  C-^  l) 

where  rii(<)  6  /2"  is  the  resultant  disturbance  vector  given  by 

ni(t)  =  nx{Aq{t),AT,^) 

=  Jv'  •  +  J'n  ■  AT{i.).  (25) 

For  some  positive  constants  £1  and  £2  defined  in  (10c),  let  the  constant  N  be  defined  as 
follows 

N  =  max{||  nx{Aq{t),AT{t),^)  ||;<7(0  €  B{ex\qd{t))  ^  qit)  E  8(62;  qdit))}  (26) 

where  the  matrix  norm  is  defined  as  the  induced  Euclidean  norm,  and  for  a  positive  number 
/>  >  0  and  a  vector  u  ^  W'  the  boundary  set  defined  by  as 

F(p;i/)  =  {le  6  /?";||  u;  -  ||<  p}.  (27) 

In  Equation  (25),  the  resultant  disturbances  are  mainly  dependent  on  the  acceleration  infor¬ 
mation  error  and  the  control  input  computation  delay  error  other  than  system  uncertainties 
or  modeling  errors  of  robot  manipulators.  The  disturbance  observer  can  compensate  mod¬ 
eling  errors  of  the  inertia  matrix  besides  SGD.  Thus  the  SMC  design  is  independent  of 
maximum  bound  of  modeling  errors  in  the  parameter  space  but  dependent  on  only  the  ac- 
( deration  information  error  and  the  control  time  delay  due  to  the  digital  implementation. 
Since  a  control  input  is  continuous  and  accurate  acceleration  can  be  obtained,  the  maximum 
of  resultant  disturbances,  N  becomes  very  small,  thus  we  can  design  a  new  SMC  without 
chattering. 
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Now,  a  stability  property  of  the  system  (5)  with  control  laws  (18)  -(22)  will  be  stated  in 
the  next  theorem. 

Theorem  1  :  Consider  the  robot  system  given  by  Equations  (18)-(22).  Assume  that  fc' 
some  positive  7,  ||  5(^0)  ||<  7  and  ||  x{to)  ||<  7/K  are  satisfied  at  the  initial  time  t  =  to,  and 
if  the  gain  k^2  satisfies 

k^2>N-k^i-6  (28) 


for  given  k^^i  and  8,  then  the  global  control  system  is  uniformly  bounded  (i.e.  the  solution 
A''  is  uniformly  bounded  at  origin  in  state  space)  for  all  t  >  to  except  ||  s(f)  ||>  rji  where  r/i 
is  defined  by 

-  [6  +  ^ ^ 

r/1  =  +  f3f  -  Qi,  ai - ^ - ,  pi  =  -j—.  (2!J) 

Proof 

If  we  take  V(f)  =  l/2s^{t)s{t)  a.s  a  Lyapnuov  candidate  and  differentiate  with  respect  ti> 
time,  it  follows 

=  s^{t)  ■  s 

=  ■  ni{t)  -  s^(0  •  {k^i  •  s{t)  +  k^2  ■  i'^O) 

[f  we  use  the  matrix  inequality 


||<||  X  II  •  |j  ?/  II  ■  II  A  II 


it  follows 


<11  S{t)  ll-llnao  II  ^(011' 


^it)  IP 


(••11) 


dt  -X.  M  -V-.  -I  -X* 

So  that  II  ni(f)  ||<  N  is  satisfied  from  the  definite  of  N,  and  we  can  rewrite  Equation  (-'ll) 


as 


<11  ^(0 II  (iV  -  II  5(1)  II  -t,2  ■ 


dt 


s{t)  II  -k^i 


II  4t)  II  +<5- 

s{t)  11^  -f2  •  ai-  II  s{t)  II  —Pi} 


II  s{t)  II  -|-(J 

at  t  =  ti.  If  the  gain  ky^i,  and  ^^2  satisfy  the  condition  (27), 


(32) 


dVjt) 

dt 


<  0 


(33) 


at  t  >  ti  until  11  s{t)  1|>  yjal  -t-  /??  —  Qi  which  completes  the  proof  of  Theorem  1. 

Theorem  1  guarantees  the  uniform  bounded  stability  of  the  proposed  continuous  SMC 
(18)-(22)  for  robot  manipulators.  This  controller  is  shown  in  Fig.  1,  its  structure  is  simple. 
The  sampling  time  can  be  as  small  as  possible  so  that  the  acceleration  information  is  dcctirat(.' 
enough,  therefore,  the  maximum  values  N  is  very  small.  If  a  smaller  6  is  used  in  control 
algorithm  (22),  then  the  lower  bound  of  rji  is  decreased.  The  lower  bound  rji  can  be  decreased 
by  increcise  of  k^i  for  given  6  and  N  so  that  r;,  is  sufficiently  small  in  comparison  to  7.  If 
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the  initial  position  error  is  small  which  is  reasonable  in  case  of  the  known  initial  state  of 
robot  manipulators,  it  is  evident  that  s(t)  remains  closer  to  the  sliding  surface  s{t)  =  0  and 
the  trajectory  error  is  also  small,  hence  we  may  assume  smcill  7  for  the  desired  specifications 
of  the  sm£ill  trajectory  errors.  The  feedback  control  (18)-(22)  designed  by  Theorem  1  and 
Lemma  1  maintains  the  stability  of  the  system  with  the  prescribed  performance.  To  design 
the  proposed  sliding  mode  controller,  firstly  the  desired  sliding  surface  defining  which  defines 
the  desired  error  dynamics  is  chosen,  i.e.,  the  gains,  and  Kp  are  specified.  Then  the  gains, 
and  k^2-i  in  Equation  (22)  are  selected  based  in  Theorem  1.  The  design  does  not  need 
the  information  of  maximum  bound  of  system  parcimeter  variations  or  uncertainties  because 
of  the  efficient  on-line  compensation.  The  proposed  SMC  can  be  designed  to  be  accurate 
and  robust  against  to  parameter  uncertainties. 


qd(t) 

qd(l) 


Fig.  1.  The  block  diagram  of  the  proposed  SMC. 


0(1)  = 


s(t) 

lls(i)li  +  8 
s(t) 

ls(i)l  +  5 


for  Theorem  1 

for  Theorem  2 


4  Modified  Sliding  Control  for  Robot  Manipulators 

In  this  section,  we  present  a  modified  control  algorithm  from  the  control  law  of  the  previous 
section  so  that  the  efficient  computation  is  possible  in  the  case  of  the  digital  implementation. 
The  modification  replaces  in  Equation  (22)  by 

T^{t)  =  -Js  ■  (k^i  ■  3{t)  +  k^2  •  <T6(f))  (of) 
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where  the  i  -  th  element  of  vector  =  [cri(t),  <72(t),  •  •  ■ ,  is  defined  by 


(35) 


I  •Sf(0  I 

If  we  apply  the  control  given  by  Equations  (18)-(21)  and  (34)  to  the  system  (1),  then 

s{t)  =  ni{t)  -  ky.2  ■  s{t)  -  ky.1  ■  ab{t)  (36) 

where  the  disturbance  vector  ni(t)  is  also  given  by  Equation  (25).  Let’s  denote  its  i  —  th 
element  by  For  some  positive  numbers  7,  e  and  >  0,  let  the  constant  N^o  be  defined 

as  follows 

=  sup  (X]  II  nu{Aq{t),AT{t),$  ||;(7(f)  €  E{£i;qci{t)  A  q{t)  e  E{£2:qit))}  (37) 

t,z,w,y,  . 


where  the  neighborhood  set  is  defined  as 

E{p;  i/)  =  {w  €  iZ";  ||  lioo<  p} 


(38) 


for  any  scalar  p  and  vector  1/  €  fC',  where  ||  •  ||oo  is  the  infinity  norm  of  the  vector  which  is 
defined  as  the  maximum  absolute  value  of  its  components. 

Theorem  2  :  Consider  the  robot  system  (5)  with  controls  given  by  (18)-(21)  and  (31). 
.Assume  that  for  some  7  >  0,  H  s(fo)  ||oo>  7  and  ||  x{to)  ||c«d<  7/«  is  satisfied  at  the  initial 
time  tg.  If  the  gain  L\2  satisfies 

kx2  <  ^00  -  6  ■  ky.1  (3!)) 

for  given  ky^i  and  S,  then  the  global  control  system  is  uniformly  bounded  (i.e.  the  solution 
.V  is  uniformly  bounded  at  origin  in  state  space)  for  all  t  >  tg  until  ||  s{t)  ||<  772  where  r/o  is 
defined  by 

-  ^  ^  I  I  ^ 

72  =  \Jotl  +  /?2  -  ^2,  Ol  =  - ^  ■  (10) 


Proof 

If  we  take  V[t)  =  l/2s^(<)s(f)  as  a  Lyapnuov  candidate  and  differentiate  it  with  respect 
to  time,  it  follows 


dV{t) 


dt 


=  •  s(0 


=  {s^  •  n,(f)  -  s^(<)  •  {ky^i  ■  s{i)  +  k^2  ■  s[t)  ■  cr(,(t))}. 
We  can  rewrite  the  above  equation  as 


dV(t) 

dt 


=  {^s,(0  ■  n,.(0  ^.(0^  -  11^x2 


So  it  follows 

dV{t) 

dt 


<  I  •3.(0  I  •  I  ”1.(0  5.(0^  -  H  ^x2 


■3j(0  I 
sAtr 


}• 


(ir 


(12) 


I  ^3,(0  I  +.5 


}■ 


>13) 


9 


From  the  definition  of  Ncx>  and  norm  inequality  of  ||  s(f)  ||>1|  s{t)  ||oo,  we  can  rewrite  the 
equation  (43)  as 

<11  S{t)  lU  {iVoo  -  II  s{t)  Iloo  -k^2  ■  II 


c^r 


s{t) 


II  s{t)  lU  +<5 

If  the  gain  k^i,  and  and  ^^2  satisfy  the  condition  (39), 


•*(0  IlL  +2  •  012-  II  s{t)  II2  -f32}. 


(44) 


dVjt) 

dt 


<  0 


(4.5) 


Cit  t  ^  until  II  s(i)  ||>  —  0(2,  the  proof  of  Theorem  2  is  completed. 

The  desired  performance  of  the  modified  controller  for  robot  manipulators  can  be  ob¬ 
tained  by  using  Theorem  2.  The  feedback  sliding  mode  control  (18)-(21)  and  (34)  maintain 
the  uniformly  bounded  stability  of  the  system  as  far  as  the  condition  (39)  is  satisfied,  but 
Theorem  2  is  only  a  sufficient  condition  for  stability.  Under  the  same  conditions,  defined 
in  Equation  (37)  is  always  larger  than  N  defined  in  Equation  (26).  Thus  as  shown  in  Equa¬ 
tions  (39)  and  (40),  for  the  same  trajectory  error  the  gains  of  the  modified  controller  are 
designed  to  be  larger  than  those  of  the  SMC  defined  in  section  3.  Under  the  same  conditions 
for  the  both  controls,  the  error  of  the  modified  version  be  smaller  than  that  of  the  original 
SMC  because  of  ||  s{t)  ||>|  3i{t)  |. 


5  Numerical  Simulation 


The  numerical  simulations  are  performed  for  the  purpose  that  is  to  show  accurate  and 
robust  trajectory  tracking  property  of  the  proposed  algorithms  compared  with  the  CT.M. 
The  SCARA-type  two  degree-of- freedom  manipulator  dynamics  [1]  used  for  simulation  are 
as  follows 


T\ 

T2 


jmi  4-  jm2  +  m2C2  3^2  "I"  \'^2C2 

qi 

\m2  -k  \Tn2C2  \m2 

.  ^2  . 

+  l2- 


—  |m252<72  ~  '^2*5291  ?2 

\m2s2qi 


+  /  ' 


^rriigCi  +  +  ^29^1 

\m2gC12 


(46) 


where  C,,  5,  and  Cq  imply  cos(gi),  3in{qi)  and  cos{qi  -f-  qr^),  respectively.  The  parameters 
are  mi  =  m2  =  0.782[%],  /  =  0.23[m]  and  g  =  9.8[m/sec^].  The  reasonable  modeling  error 
of  robot  dynamics  is  produced  for  each  ma.ss  and  length.  Moreover,  unknown  payload  is 
reflected  to  m2.  To  test  the  robustness  to  modeling  uncertainties  and  unknown  payload, 
simulation  studies  are  carried  for  the  CTM  and  two  proposed  SMCs  under  three  case  condi¬ 
tions,  i.e.,(i):no  modeling  error  and  no  unknown  payload,  (ii):10  [%]  modeling  error  and  0.5 
[kg]  unknown  payload,  and  (iii);  10  [%]  modeling  error  and  1  [kg]  unknown  payload.  The 
desired  trajectory  for  each  link  wais 


+ 


(<7/-<7.)  ,  {qj  -  qi)  ■  sin{Trt)  ^  180  , 

— 2 - < - - - }.—  [degree] 


(47) 
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where  the  initial  position  9,  =  [0.0  0.0]^  [degree],  and  the  final  position  qj  =  [60.  -  60.]^ 
[degree].  We  take  the  execution  time  of  the  trajectory  tracking  as  2  seconds.  The  necessary 
time  for  the  computation  of  model  (2)  was  selected  to  be  2  [msec],  and  this  computation  time 
is  considered  as  a  time  delay  of  input  for  the  simulation  of  the  CTM.  The  gain  kp  =  400  •  /. 
ku  =  40  •  /  were  used  in  the  sliding  surface  (6)  so  that  the  ideal  sliding  dynamics  have 
its  eigenvalues  to  be  assigned  to  -20.  The  corresponding  constants  defined  in  (9)  ,  K  and 
K  become  14.7  and  10.  Using  the  results  of  Lemma  1,  the  trajectory  error,  .Yi,  and  its 
derivative,  X2,  are  bounded  as  Ci  =  1.47  ■  7  and  €2  =  592  •  7  for  a  given  7  of  the  constraint  to 
the  sliding  surface.  The  smaller  bound  of  the  sliding  surface,  7,  the  smaller  trajectory  error 
and  its  derivative.  For  0.01  [degree]  maximum  tracking  error,  7  was  selected  to  be  0.006S. 

Now,  in  order  to  obtain  the  above  prescribed  tracking  performance  ,  the  propo.sed  SM(  's 
are  designed  to  make  the  sliding  surface  be  bounded  by  7.  Assuming  the  initial  position 
of  the  robot  manipulator  is  known,  the  conditions,  ||.Y(to)||  <  7/«:  and  ||s(to)||  <  *  arc 
satisfied  at  initial  time  t  =  to-  The  controller  gains  k^^i  and  ky.2  were  selected  to  be  50  and 
10  which  satisfy  the  conditions  (39)  and  (45)  for  given  8  =  0.05  and  N  or  so  that  r/i 
and  772  are  sufficiently  small  with  respect  to  the  selected  7  by  Theorem  1  and  Theorem  2. 
For  comparison  of  the  three  case  conditions,  the  N  and  Noo  are  summarized  in  Table  1.  .\s 
the  payload  increases,  the  value  N  and  Noo  increased.  As  a  results  the  lower  bound  of  k\2 
becomes  larger  such  that  the  stability  of  the  system  is  maintained.  The  simulation  block 
diagram  of  proposed  algorithms  are  shown  in  Fig.  1.  The  proposed  trajectory  tracking 
controller  is  consisted  of  the  modified  equivalent  control,  feqit)  (21),  feedback  of  the  sliding 
surface,  T^{t)  (22)  and  disturbance  observer,  Tc{t)  (19),  terms. 

The  gains  A'y  =  40  •  /,  Kp  =  400  •  I  were  used  in  the  CTM  algorithm.  The  simulation 
results  of  the  CTM  are  shown  in  Fig.  2.  The  tracking  errors  for  the  CTM  are  shown  in 
Fig.  2-1  (a)  for  link  one  and  (b)  for  link  two,  and  its  corresponding  control  inputs  are  shown 
in  Fig.  2-2  (a)  for  link  one  and  in  Fig.  2-2  (b)  for  link  two.  As  the  payload  increases, 
the  root- mean-square  (rms)  error  is  rapidly  increased  from  0.0073  [degree]  to  3.6801  [degree] 
using  the  CTM  algorithm.  Thus  the  CTM  need  to  be  added  the  robust  algorithm  against 
model  uncertainties  and  unknown  payloads. 

The  results  of  the  proposed  SMC  simulations  are  shown  in  Fig.  3  for  Theorem  1  and 
in  Fig.  4  for  Theorem  2.  The  tracking  errors  of  the  proposed  SMC  are  presented  in  Fig. 
3-1  and  Fig.  4-1.  The  tracking  errors  by  Theorem  2  are  smaller  than  those  by  Theorem 
1  because  the  norm  of  the  sliding  surface  in  (22)  is  greater  than  the  infinity  norm  of  the 
sliding  surface  in  (35).  The  correspondent  continuous  control  inputs  and  sliding  surfaces  are 
depicted  in  Fig.  3-2  and  Fig.  3-3  for  Theorem  1  and  in  Fig.  4-2  and  Fig.  4-3  for  Theorem 
2.  The  simulation  results  in  view  of  maximum  and  rms  values  are  summarized  in  Table  1. 
Using  the  proposed  algorithm,  the  resultant  error  maintains  nearly  unchanged  Especially, 
note  that  using  the  proposed  algorithm  the  tracking  error  is  not  increased  so  much  a.s  the 
increase  of  the  palyload  error,  for  the  case  Theorem  1  and  Theorem  2. 

From  the  simulations, we  have  found  that  the  proposed  two  algorithms  provides  a  better 
performance  than  the  CTM  algorithm  with  respect  to  the  tracking  errors.  Moreover,  the 
tracking  errors  of  the  proposed  algorithms  satisfy  the  prescribed  specifications  as  shown  in 
Fig.  3-1  and  Fig.  4-1  subject  to  parameter  and  payload  uncertainties. 
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Table  I.  Numerical  comparison  between  CTM  and  proposed  SMC 


Method 

CASE 

Error 

(max) 

[degree] 

Error 

(rms) 

[degree] 

Control 

(max) 

[Nm] 

Control 

(rms) 

[Nm] 

Surface 

(max) 

Surface 

(rms) 

N 

CTM 

1 

0.0066 

0.0073 

3.6438 

3.7569 

2 

3.6375 

3.6801 

5.7958 

5.8300 

3 

6.3872 

6.4713 

8.0360 

8.0453 

SMC 

THM  1 

1 

0.0004 

0.0005 

3.6438 

3.7569 

0.0015 

0.0016 

0.3281 

2 

0.0008 

0.0010 

5.3968 

5.6014 

0.0031 

0.0034 

0.4862 

3 

0.0011 

0.0016 

7  Am 

7.4476 

0.0041 

0.0044 

0.6462 

THM2 

1 

0.0002 

0.0002 

3.6439 

3.7569 

'  0.0008 

0.0008 

0.4681 

2 

0.0004 

0.0005 

5.3968 

5.6014 

0.0015 

0.0016 

0.6911 

3 

0.0005 

0.0006 

7.1498 

IMlb 

0.0020 

0.0021 

0.9153 
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TRACKING  ERROR  2  (degree)  TRACKING  ERROR  1  [degree] 


k 


10  -r 

9  -i 


-4  '  I  T  t  t  r  1  1  I  I  \  . .  . f  ■  '-I 

0  0.2  0  4  0.6  0.8  1  1  2  1  4  1,6  1  8  2 


TIME  [sec] 

(b)  link  two 

Fig.  2-1  Tracking  errors  for  CTM  algorithm 

(i)  case  1  :  no  model  error  and  no  unknown  payload 

(ii)  case  2  :  10  [%]  model  error  and  0.5  [kg]  unknown  payload 

(iii)  case  3  :  10  [%]  model  error  and  1.0  [kg]  unknown  payload 
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CONTROL  2  [NmJ  CONTROL  1  [Nm] 


10 


(i)  :  CASE  1 

(ii)  ;  CASE  2 

(iii)  :  CASE  3 


0  2  0  4  0  6  0.8  1  1.2  1.4  1  6  1  8 


TIME  [sec] 
(a)  link  one 


TIME  [sec] 
(b)  link  two 


Fig.  2-2  Corresponding  control  inputs  by  CTM  algorithm 
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TRACKING  ERROR  2  (degree]  TRACKING  ERROR  1  [degree] 


TIME  [sec] 
(a)  lini:  one 


TIME  [sec] 
(b)  link  two 


Fig.  3-1  Tracking  error  for  Theorem  1  of  proposed  SMC 

(i)  case  1  :  no  model  error  and  no  unknown  payload 

(ii)  case  2  :  10  {%]  model  error  and  0.5  [kg]  unknown  payload 

(iii)  case  3  :  10  [%]  model  error  and  1.0  [kg]  unknown  payload 
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CONTROL  2  (Nm)  CONTROL  1  (Nm) 


i 


TIME  [sec] 
(a)  link  one 


(i)  :  CASE  1 

(ii)  :  CASE  2 

(iii)  :  CASE  3 
(iii) 

- _ 01 _ _ _ 

0  02  04  06  06  1  1  2  1  4  1  6  1  8  2 


TIME  [sec] 

(b)  link  two 

Fig.  3-2  Corresponding  control  inputs  by  Theorem  1  of  proposed  SMC 
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SLIDING  SURFACE  2  SLIDING  SURFACE  1 


0.005 


TIME  [sec] 

(b)  link  two 

Fig.  3-3  Corresponding  sliding  surfaces  for  Theorem  1  of  proposed  SMC 
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TRACKING  ERROR  2  (degree)  TRACKING  ERROR  1  (degree) 


(a)  link  one 


i 


(b)  link  two 

Fig,  4- 1  Tracking  errors  for  Theorem  2  of  proposed  SMC 

(i)  case  1  :  no  model  error  and  no  unknown  payload 

(ii)  case  2  :  10  [%]  model  error  and  0.5  (kg)  unknown  payload 

(iii)  case  3  :  10  [%]  model  enror  and  1.0  (kg)  unknown  payload 
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CONTROL  2  (Nm)  CONTROL  1  [Nm] 


0  02  04  06  08  1  12  14  16  18  2 

TIME  [sec] 

(a)  link  one 


t 


0  02  04  06  08  1  1  2  1  4  1  6  1  8  2 

TIME  [sec] 

(b)  link  two 

Fig.  4-2  Corresponding  control  inputs  by  Theorem  2  of  propoed  SMC 
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SLIDING 


0  004  -1 


(i)  ;  CASE  1 

(ii)  :  CASE  2 

(iii)  :  CASE  3 


0  003  - 


LU 


■0  00*  ^ 


-0  006  . . .  I  '  1  '  '  !  !  t  ■  1  ! - r— 

0  02  04  06  08  1  12  14  16  18 

TIME  [sec] 

(a)  link  one 
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Fig.  4-3  Corre.sponding  sliding  surface  for  Theorem  2  of  proposed  SMC 
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6  Conclusions 


In  this  paper,  continuous  sliding  control  algorithms  have  been  proposed  for  accurate  and 
robust  tracking  control  of  robot  manipulator  and  their  stability  properties  are  analyzed  in 
detail.  As  a  preliminary  work,  the  relationship  between  the  maximum  bound  of  the  tracking 
error  and  that  of  the  sliding  surface  is  presented  in  Lemma  1.  The  continuous  sliding  mode 
control  algorithm  with  disturbance  observer  is  proposed  and  its  uniform  bounded  stability 
property  is  proven  in  Theorem  1.  To  improve  efficiency,  a  modification  to  the  proposed 
algorithm  has  been  preceded  and  its  stability  property  is  presented  in  Theorem  2.  Unchn- 
these  control  algorithms,  the  tracking  errors  can  be  reduced  to  the  prescribed  range  without 
chattering  problems  based  on  the  above  stability  analysis. 

The  simulation  results  have  shown  that  the  proposed  algorithms  are  very  efficient  for  the 
trajectory  tracking  and  robust  to  the  modeling  inaccuracy  and  unknown  payloads. 
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